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ABSTRACT 

Context. The Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory (ISO) observed the four large main-belt 
asteroids (1) Ceres, (2) Pallas, (4) Vesta, and (10) Hygiea multiple times. The photometric and spectroscopic data cover the wavelength 
range between 43 and 197 /jm, and are a unique dataset for future investigations and detailed characterisations of these bodies. 
Aims. The standard ISO archive products, produced through the last post-mission LWS pipeline, were still affected by instrument 
artefacts. Our goal was to provide the best possible data products to exploit the full scientific potential of these observations. 
Methods. For all asteroid observations we analysed in detail the dark current, the calibration reference flashes, the space environment 
effects (glitches), memory effects, tracking influences, and various other sources of uncertainty. We performed a refined reduction 
of all measurements, corrected for the various effects, and re-calibrated the data. We outline the data reduction process and give an 
overview of the available data and the quality of the observations. We apply a thermophysical model to the flux measurements to 
derive far-IR based diameter and albedo values of the asteroids. The measured thermal rotational lightcurve of (4) Vesta is compared 
to model predictions. 

Results. The catalogue of LWS (Long Wavelength Spectrometer) observations of asteroids contains 57 manually reduced datasets, 
including seven non-standard observations, which as such did not have final pipeline products available before. In total, the archive 
now contains 1 1 spectral scans and 46 fixed grating measurements with a simultaneous observation at 10 key wavelengths distributed 
over the full LWS range. The new data products are now accessible via the ISO data archive as highly processed data products 
(HPDP)**. 

Conclusions. The quality of the data products was checked against state-of-the-art thermophysical model predictions and an excellent 
agreement was found. The absolute photometric accuracy is better than 10 %. The calibrated spectra will serve as source for future 
mineralogical studies of dwarf planets and dwarf planet candidates. 

Key words. Minor planets, asteroids - Radiation mechanisms: Thermal - Infrared: Solar system 



q ; 1 ■ Introduction ISO Handbook, Volume III dGry et alJl2003h . In total, the ISO 
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> , The Infrared Space Observatory (ISO; [Kessler et al. | [19961) of four different asteroids coveri almost 18 hours of satellite 

• . observed between 1995 and 1998 more than 40 asteroids in ^ However, the scientific return of these programmes was so 

; great detail, including some complete spectra fro m 2 to 200^m far lhtle Most of the scientific quest ions could not be an- 
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face alteration processes, and the interpretation of taxonomic ..° „ .... .... , . ° r . . . 

. ^ 11J -.C- r • j • r jj? r future the Herschel mission will be capable of obtaining aster- 
classes through the identification of mid-infrared features of . , . ,_ , f, ,™ ? 
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. , . . „„ . j. . : — — , — Array Camera & Spectrometer) and between 200 and 670 urn 
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, „ „ , T „„ . J r with SPIRE (Spectral and Photometric Imaging Receiver), but 

grammes and from all tour ISO instruments. ., j j, 

m iit i o „,,,„ — n asteroids are currently not part of the accepted key programmes. 

The ISO Long Wavelength Spectrometer (LWS; Clegg et al. .... ... TC „ . J ,. ./ A • , . * 

I , a , , . • j • Although the ISO target list is very limited, it includes obiects 

1996) performed photometric and spectroscopic measurements c ,. , . , . ° . . . 

. — -r" r . ^, , n , m j . .. of high scientific interest: (1) Ceres and (4) Vesta are the targets 

in the wavelength range 43 to 197 um and opened, together with ,., e _ ra a. • • , ( >n . , ., . , 

mnminr ■ j r r ■ r j j i of the Dawrli spacecraf t mission, ( 1 ) Ceres is mean while consid- 

ISOPHOT, a new window for far-infrared astronomy. The details , , _J , j^nn . , ,. mlI 

. . ■ ■ , • i^i-, ered as a dwarf -planet, and (2) Pallas, (4) Vesta, and (10) Hygiea 

on the instrument and the data processing can be found in the , „ . f ... . h T , .. . , ,. , ./ 

r ° are dwarf-planet candidate^. In general, it is believed that these 
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largest bodies in the main-belt are protoplanets remaining intact 
since their formation. 

The LWS observations comprise spectral energy distribu- 
tions over the full wavelength range between 43 and 197 /mi with 
a spectral resolving power of ~ 200 and can be utilised to search 
for spectral signatures related to the asteroids' mineralogic com- 
position. An overview of the spectroscopic features from theory 
and lab measurements in t he infrared wavel ength regime up to 
25 /mi can be found in e.g. S alisbury ( 1993) and in the context 
of ISO in Pot to et al.l (120021) . Far infrared laboratory spectra of 
various crystalline and amorphous minerals are digitally avail- 
able e.g. in the Heidelberg - Je na - St.Petersburg - D atabase of 
Optical Constants (HJPDOCfl (H enning et aill 19991) . Although 
these measurements are in principle sufficient to allow the anal- 
ysis of the data presented in this paper, their application to solid 
surfaces with non-uniform grain sizes and subsurface radiation is 
not straightforward, rendering a mineralogic interpretation diffi- 
cult. 

Apart from that, the LWS observations can give valuable in- 
put to models describing the thermal emission of minor bodies in 
our own or other solar systems, and are key ingredients to estab- 
lish a steroids as calibration objects in the middle and far infrared 
(e.g. lMiiller & Lagerroslll998ll2002h . 

An effort was made to produce a homogeneously reduced set 
of these observations, eventually leading to a catalogue which is 
available through the IDA. We describe the contents of the cata- 
logue (Section|2]), outline the data reduction process, and give an 
overview of the instrumental artefacts (Section [3]) encountered 
during catalogue compilation. Section [4] gives a brief overview 
of the catalogue structure. In Section [5] we discuss the quality 
of the final products, present diameter and albedo values derived 
from thermal IR data in the range 43-197/mi for the observed 
asteroids, and analyse the thermal lightcurve of (4) Vesta. 

2. Observations 

2.1. Overview 

Most of the LWS asteroid observations were carried out in fixed 
grating mode, using the astronomical observing template (AOT) 
L02. This produced ten photometric measurements - one in 
each detector - at 46.2 /mi, 56.2 //m, 66. 1 //m, 75.7 /mi, 84.8 /mi, 
102.4/mi, 141.8/mi, 160.6//m, and 178.0/mi (corresponding to 
detectors SW1 to SW5 and LW1 to LW5). 

Only four full grating scans using the template L01 were per- 
formed, covering the whole spectral range from 43 to 197 /mi 
with typical SNRs of 20-30 for (10) Hygiea and 150-200 for 
(1) Ceres. The spectrum is composed of ten sub-spectra with the 
sub-spectra being generated by the grating scanning over the ten 
LWS detectors simultaneously. Signal-to-noise was built up by 
taking more than one scan. 

Additionally there exist 25 observations carried out in the 
non-standard engineering mode L99, most of them similar to 
L01 grating scans (SNR * 150 for (1) Ceres and ^100 for (4) 
Vesta). For the catalogue, only L99 on-source grating scans were 
considered, and no background measurements are included. The 
non-standard observations are characterised by experimental set- 
tings of the detector parameters (bias voltages and heater cur- 
rents), an unusually high spectral sampling (13 samples per res- 
olution element instead of four, as it was the standard setting 
later), and longer integration times per grating position. 

From the available L99 observations we selected seven 
datasets for inclusion in this catalogue. The spectral analysis 

4 http : //www . mpia-hd . mpg . de/HJPDOC/ 



of this data still remains difficult, since the relative spectral re- 
sponse function (RSRF) of the detectors used in the calibration 
process is based on observations performed in the standard L01 
mode. The remaining L99 data sets did not include sufficient cal- 
ibration data for reliable dark current determination. 

All standard observations are listed in Tbl. [T] while the non- 
standard observations are summarised in Tbl. [2] Contamination 
by 158 /mi CII background emission is indicated where applica- 
ble. 

A comprehensive overview of the instrument design, observ- 
ing modes , and detector p roperties can be found in the LWS 
handbook dGrv et alj|2003l) . 

2.2. Observation elements 

Each observation consists not only of flux measurements of the 
target itself, but also includes calibration data, typically split in 
two blocks at begin and end. By placing the Fabry-Perdt spec- 
trometer into the beam and intentionally misaligning the etalons, 
light from the source did not reach the detectors, eventually re- 
sulting in a dark current measurement. Long observations, and 
especially the engineering mode L99 data, contain intermittent 
additional dark current measurements. 

Since the responsivity of the detectors is not constant, but 
usually drifting upwards over time due to impacts of charged 
particles, the actual responsivity of each detector had to be deter- 
mined for each observation separately. This was done by flash- 
ing the detectors with internal illuminators while light from the 
source was again blocked with the Fabry-Perdt. This allowed to 
quantify responsivity changes during the course of a single ob- 
servation, and to apply the so-called 'responsivity drift correc- 
tion' . 

Absolute responsivity calibration was performed by compar- 
ing the responsivities measured at the time of the observation to 
the ones measured during observations of Uranus, the primary 
LWS calibrator. 

For the target source flux measurements, the Fabry-Perdt was 
removed from the beam. Depending on the observation template 
and requested wavelength coverage, the grating remained either 
at a fixed angle or was moved to scan the desired spectral range. 
This could be done in unidirectional mode, where the grating an- 
gle constantly in- or decreased and was reset to its start position 
before each scan repetition. Alternatively, the observation could 
be performed in bidirectional mode, with the grating angle in- 
creasing and then decreasing with the same rate during the mea- 
surement. This is referred to as up- and down-scan, where 'up' 
and 'down' does not correspond to an increase or decrease of the 
observed wavelength, but of the grating position encoder value. 
In fact, during an upscan the observing wavelength decreased, 
and increased during downscans. 

3. Data reduction 

3. 1 . Instrumental artefacts 

The most prominent artefacts, common to all ISO instruments, 
are 'glitches', sudden increases of the measured flux due to 
charged particles hitting a detector or parts of the readout elec- 
tronics. The initial sudden signal increase is sometimes followed 
by a 'glitch tail' , a long-lived exponential decay. Though glitches 
can be recognised and filtered out automatically to some extent, 
glitch remnants are present in all observations throughout the 
archive, regardless of instrument or observing mode. Glitches do 
not only occur during object flux measurements, but also during 
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the calibration, i.e. dark current and responsivity measurements, 
and adversely affect the accuracy and reliability of the actual 
observation. Special attention was therefore put on the manual 
removal of glitch remnants in each observation contained in the 
presented catalogue. 

Other frequently observed instrumental artefacts are fringes 
in the spectra in the case of extended sources or off-axis ob- 
servations, spurious spectral features due to uncertainties in the 
spectral response function used in the calibration process, and 
memory effects, i.e. changes of detector response depending on 
the illumination history. 

Observations towards the end of the satellite mission were 
affected by so-called 'warm-up effects', visible as broad spectral 
features in some detectors. This was presumably caused by peri- 
odic breaks of the liquid helium film in the vicinity of the LWS 
strap location as the satellite's helium tank came close to ex- 
haustion. From the observations contained in this catalogue, this 
affects one L01 type observation of (1) Ceres taken in December 
1997, several months before ISO's liquid helium depletion in 
April 1998. 

A thorough description of caveats and unexpected effects 
concerning LWS data can be found in Chapter 6 of the LWS 
handbook dGrv et al.ll2003l) . All datasets included in this cata- 
logue were checked for the presence of the artefacts described 
above and data quality was quantified via a set of well-defined 
flags. The criteria for setting these flags and the results of this 
quality check can be foun d in the technical docume ntation ac- 
companying the catalogue (iHormuth & Miiller 2006). 

3.2. Data reduction steps 

Data reduction was performed with the available standard reduc- 
tion packages LIA (LWS Interactive Analysis) Version 10.2 and 
ISAP (ISO Spectral Analysis Package) Version 2.2. Additional 
own IDL-procedures were used for manual editing of the data at 
the SPD (Standard Processed Data) level, and the generation of 
the final FITS and ASCII files. The SPD data is preprocessed in 
the sense that raw detector readouts have been converted into 
photocurrents, saturated readouts flagged, and strong glitches 
automatically identified and removed. Housekeeping data not 
necessary for further data reduction has been removed, but dark 
current and responsivity measurements have been preserved. In 
the following we briefly outline the order of reduction steps and 
give the names of the corresponding LIA procedures where ap- 
plicable. 

In a first step, the pointing coordinates of each observation 
were checked against the ephemeris of the observed asteroid to 
detect pointing problems. The maximum difference between as- 
teroid position and telescope pointing was found to be 2 '.'5, much 
smaller than the LWS beam profile (see iLlovdl l2003h . The re- 
duction process started at SPD level with a visual inspection 
of the data in the time domain. Photometric L02 observations 
were manually de-glitched at this stage, capturing minor glitches 
not automatically flagged during generation of the SPD data. 
The available pipeline data of grating scans was additionally in- 
spected with ISAP to get a quick overview of the available scans 
and scan directions. 

Next, dark current determination and subtraction were per- 
formed using the LIA routine IAJ3ARK, followed by the abso- 
lute responsivity correction with IA_ABSCORR. These interactive 
routines allow visual inspection and manual editing of the dark 
current measurements and the internal calibration source obser- 
vations. Deglitching of this data and the removal of measure- 
ments obviously affected by memory effects is critical, since any 



artefacts in the calibration data will affect the accuracy and relia- 
bility of the whole observation. The LIA routine SH0RT_AAL was 
used to create ISAP-compatible output spectra in the LWS Auto- 
Analysis format, and to convert to W/cm 2 /yum, the flux unit used 
for all LWS observations. This routine automatically strips all 
calibration data and leaves only the basic information needed for 
further scientific analysis, e.g. wavelength, flux, detector num- 
ber, and scan direction. 

In the case of photometric L02 observations, the data was 
now averaged per detector to create the final FITS product. 

In the case of grating scans a first round of manual de- 
glitching was performed, probably the most time consuming step 
in the reduction process. Relative responsivity correction was 
achieved by comparing the mean flux in each detector and scan 
to the mean flux level of all scans together. The derived gains 
were then used to bring all scans to the same flux level - sepa- 
rately for each detector and scan direction. 

After fine-zapping of remaining glitches, the data was aver- 
aged per detector with a bin width of typically 0.06 //m, corre- 
sponding to an oversampling of four with respect to the spectral 
resolution. In the case of bidirectional scans the data was now av- 
eraged over the scan directions. Especially at longer wavelengths 
some spectra clearly showed fringing, probably caused by back- 
ground emission. The affected detectors were defringed within 
ISAP, using only the averaged spectrum for fringe detection and 
defringing to get the best possible signal-to-noise for the fringe 
fitting (in the case of bidirectional data). 



4. Data products and access to the catalogue 

The ISO Data Archive (IDAfl contains spectroscopic, imaging, 
photometric, and polarimetric measurements of more than 40 
different asteroids at wavelengths between 2 and 240/mi. A ver- 
satile Java based web-interface allows to extract observational 
data in various formats, ranging from unprocessed raw data, over 
standard pipeline processed measurements up to manually repro- 
cessed high quality data products. 

All data files associated to this catalogue can be retrieved 
by searching for the specific TDT@ number of a measurement 
as given in Tbl.[T]and|2j by searching for all LWS observations 
of a given asteroid, or by searching for LWS asteroid observa- 
tions in general. The highly processed data products, including 
overview plots, quality flags, and catalogue documentation can 
be obtained via the HPDP button left to each entry in the search 

result window of the archive i nterface. 

The catalogue documentation (IHormuth & Miiller 2006), giving 
a thorough description of the data reduction steps, is also sepa- 
rately available from the list of all HPDP data setfl 

In Fig. [T]and Fig. [2] we show logarithmic example plots for 
all observing modes. Different colours have been chosen to al- 
low distinguishment of overlapping measurements by different 
detectors. The error bars in Fig. Q] include measurement and 
calibration errors, incuding a 5% uncertainty attributed to the 
Uranus models used for absolute flux calibration. The slight mis- 
matches between measurements from different detectors at same 
wavelengths in Fig.|2]are caused by memory effects of individual 
detectors. 

The data sets are available both as FITS tables, compati- 
ble with the ISAP software package, and as plain ASCII ta- 



5 Accessible via http : //iso . esac . esa. int 

6 "Target Dedicated Time", a unique identifier for ISO observations 

7 http : //iso . esac . esa . int/ida/hpdp . php 
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Fig. 1. Examples for L02 observations. Left: fixed grating ob- 
servation of (2) Pallas. Right: short grating scan observation of 
(lO)Hygiea. 



Fig. 2. Examples for full wavelength range scans. Left: standard 
L01 observation of (1) Ceres. Right: Engineering mode L99 ob- 
servation of (4) Vesta. 



bles. Please refer to the catalogue documentation for information 
about individual FITS header key words and table columns. 



5. Discussion 

5.1. Quality of the catalogue products 

In order to study the reliability of the photometric L02 data in 
the catalogue, we compared the measured flux with predictions 
calculated by m eans of the thermophysical model (TPM; e.g. 
lLagerrosf l998. and references therein). For Tbl.[3]we calculated 
model fluxes for the specified targets and epochs and analysed 
the observation-to-model ratios. Measurements affected by de- 
tector warm-up artefacts or abnormally high dark current were 
not included in this quality check. 

The general agreement between the photometric observa- 
tions (L02) and models (Tbl. [3j is better than 5 % for Ceres, 
Pallas, and Vesta and better than 30 % for Hygiea. The higher 
uncertainties are due to dark current effects at the low flux lev- 
els. The discrepancy for Hygiea is not surprising. Due to the 
lack of suitable good quality thermal data the thermophysical 
model parameters are also very p oor, and Hygiea is consid ered 
as a low quality far-IR calibrator (Mtil ler & Lagerroslll998l) . As 
a rule of thumb, the absolute flux error is better than 10% for 
targets brighter than about 10 Jy. 



Table 3. Ratios between LWS 02 observations and TPM predic- 
tions, weighted mean values and standard deviations per detector 
(SW1, 2, 3, 4, 5 and LW1, 2, 3) and per asteroid. The photom- 
etry of the LW4 and LW5 detectors is less reliable due to dark 
current problems. 





(1) Ceres 


(2) Pallas 


(4) Vesta 


(10) Hygiea 


A 


1 1 Obs. 


4 Obs. 


14 Obs. 


2 Obs. 


[jum] 


Obs/Mod 


Obs/Mod 


Obs/Mod 


Obs/Mod 


46.4 


1.04±0.05 


1.00+0.02 


1.02+0.04 


1.21+0.05 


56.4 


1.01±0.03 


1.06+0.02 


1.01+0.04 


1.25+0.10 


66.3 


1.04±0.02 


1.10+0.04 


1.02+0.04 


1.23+0.15 


75.9 


1.02±0.03 


1.02+0.03 


1.04+0.04 


1.24+0.12 


85.0 


1.01±0.04 


1.02+0.11 


1.00+0.07 


1.39+0.10 


102.9 


1.09±0.04 


1.13+0.04 


1.24+0.06 


1.40+0.26 


122.7 


0.96±0.05 


0.95+0.03 


1.07+0.05 


1.30+0.07 


142.2 


0.89±0.04 


0.84+0.07 


1.10+0.11 


1.33+0.04 


161.0 


0.99±0.05 


0.99+0.04 






178.4 


1.01+0.10 








Mean 


1.01±0.05 


1.03+0.08 


1.05+0.08 


1.27+0.14 



The high quality of the spectral scans can be seen in the 
smooth and well-connected curves of the archive browse images 



F. Hormuth and T. G. Miiller: Catalogue of ISO LWS observations of asteroids 



5 



(see Fig. |2J. For mineralogic studies it might still be useful to 
investigate the up- and down-scans separately to confirm that 
certain low-level, broad-band structures in the scans are reliable 
or not. Low level features extending over two or more detectors 
as well as structures close to pronounced RSRF changes have 
to be analysed carefully: the detectors are not always behaving 
in exactly the same way and LWS should be considered as 10 
individual spectrometers, one per detector. 

5.2. Diameter and albedo calculations 

Direct size measurements, even for the large main-belt asteroids, 
are still difficult. Up to now, only a small sample has been re- 
solved via imaging techniques from space (HST, see lDotto et al.l 
(2002) and references therein) and adapti ve optics technique s 
from large ground-based telescopes (e.g. Marchis et al. 2006). 
Occultation observations performed by well-organised groups 
of amateur astronomers also reveal direct sizes, expressed in 
chord maps which represent the cros s section of the asteroid 
at the time of the occultati on (e.g. iMillis & Dunhaml 11989) 
or http://www.psi.edu/pds/resource/occ.html for more recent re- 
sults). But the largest number of astero id sizes are originating 
from thermal infrared observations (e.g. Tedesco et al. 20 02albl) 
via radiometric methods i.e., an indirect method which relies on 
model assumptions. 

We used a well-establ ished thermophysical model (TPM; 
lLagerrosI 19961 [l997ll 19981) to derive absolute effective sizes and 
albedos for our 4 asteroids from the new L02 data sets. 

This model allows to use state-of-the-art shape models to- 
gether with corresponding spin axi s solutions derived from 
lightcurve inversion techniques (e.g. [Kaasalain en et al.l [2002). 
As part of the energy equation for the asteroids the amount of 
reflected sunlight i s described via the H-G magnitude system 
dBowell et al.ll 989). The rest of the solar energy is absorbed and 
re-emitted as thermal emission. Knowing the solar insolation to- 
gether with the disk-integrated thermal emission allows then to 
solve the e nergy equation for the a steroid's effective size and 
albedo (e.g. Harris & Lagerros 2002). The true illumination and 
observing geometries are taken into account in the TPM calcu- 
lations. The surface roughness is described by the r.m.s. values 
of the surface slopes (p) and the fraction of the surface covered 
by craters (f). For the surface roughness, as well as for the ther- 
mal properties to calculate the heat conduction into lower lay- 
ers of the surface, we used "default properties" (p=0.7, f =0.6, 



15 J irT 2 s- (X5 KT 1 ) as described by Miiller et al. ( 1999). The 



emissivity models were taken from Miiller & Lagerros (1998): 
a standard wavelength-dependent emissivity model for Ceres, 
Pallas, ad Hygiea and a special one for Vesta where lower emis- 
sivities in the submm/mm range were found. 
The TPM with the "default thermal properties" has been tested 
and validated extensively in the c ontext of ISO for large , 
regolith-covered main-belt-asteroids (Miiller & Lagerros 2002). 
It is worth to note here that the surface roughness and heat con- 
duction properties are less important for our far-IR analysis since 
they mainly influence the mid-IR r adiation at the Wien-part of 
the spectral energy distribution (e.g. lMullerll2002l Figs. 2&3). 
Asteroid-specific model input parameters and the corresponding 
references are listed in Tbl.|4] 

Our derived diameter and albedo values are shown in Tbl. [5] 
The diameters are effective diameters of a sphere of equal vol- 
ume. It can be considered as a "scaling factor" for shape models 
which lack the absolute size information. These values are the 
weighted mean values, based on the photometrically more reli- 
able short wavelength detectors (46.4 to 85.0/im). We also ex- 



cluded the data which are flagged as being affected by high dark 
currents or high glitch rates. In the table we give the full stan- 
dard deviation of all N radiometric solutions to account for the 
fact that different observations see the asteroid under different 
observing and illumination geometries. 



Table 5. Radiometric diameter and albedo values, derived from 
the LWS L02 measurements. 



Object 




[km] 


Pv 




N 


(1) Ceres 


959.6 


± 16.9 


0.096 ± 


0.003 


55 


(2) Pallas 


534.4 


± 14.7 


0.142 ± 


0.009 


20 


(4) Vesta 


548.5 


± 12.5 


0.317 ± 


0.015 


70 


(10) Hygiea 


469.3 


±26.5 


0.056 ± 


0.006 


20 



These diameter and albedo values are the first ones which 
are derived purely from far-IR observations. Nevertheless, the 
results agree very well with the so far published values. 

he most accurate val ues for the diameter of (1) Ceres come 
from HST observations dThomas et al.ll2005b and adaptive op- 
tics (A O) observations using the Keck II telescope dCarrv et al.l 
2008). The HST observations led to an oblate spheroid with 
axes of a=b=487.3km and c=454.7km, resulting in an effective 
diameter of 2 x (a x b x c) 1/3 =952.4km. The AO data con- 
firmed the oblate spheroid, but with slightly different values (a 
= b = 479.7+2.3 km, c = 444.4+2.1 km) and an effective diame- 
ter of 935.3 km. Our result is derived from very limited viewing 
geometries, but it is in excellent agreement with the HST re- 
sults. In this context it is interesting to note that the diameter 
derived from IR AS observations at 12, 25, 60 and lOO/vm was 
848.40+19.7km dTedesco et a l. 2002b). Also our albedo value 
agrees very well with the 0.0936 derive d via stellar occultation 
techniques (Shevchenko & Tedesco 2006). 

(2) Pallas: iTedesco et all d2002bl) obtained from IRAS 
of 498. 07+18. 8km and a geometric albedo of 
A combination of speckle, occultation and 



data values 
0.1587±0.013. 

lightcurve observations of (2) Pallas lead to an ellipsoidal shape 
with dimensions of 574 ± 10 x 526 ± 3 x 501 ± 2 km and 
a mean diameter of 533 ± 6 km dDunham et alJ 1 19901). A re- 
cent h igh-angular resolution adaptive optics study bv fCarrv et all 
d2007l) resulted in an absolute effective size of 519.5 km (semi- 
major ellipsoidal axes values of 276x256x248km, ±10km). 
IDrummond & Christoul (120081) found 494+57 km based on an- 
other set of AO observations. The various results show that there 
are still significant uncertainties in the overall shape and size 
values, but the derived effective sizes agree within the given er- 
rorbars with our radiometric diameter. Therefore, we believe that 
our albedo of 0. 142 is the most reliable value so far published for 
Pallas. 

(4) Vesta has a very complex surface w ith large albedo vari- 
ations on the surface dBinzel et al.lll997l) . Our diameter value 
is about 4% larger than the effective diameter from HST ob- 
servations (axes 289 , 280, and 229 + 5 km; D e ff = 529.2km; 
Thomas et al. ( 1997), note added in proof) and about 7% larger 
than the mean triaxial ellipsoid solution with 563 + 5 x 534 ± 
5 x 442 + 7km (2 • (q ■ b ■ c ) 1/3 = 510.3 + 5.6km) given by 
IDrum mond & Christoul (120081) . This is outside the given rms- 
scatter from our 70 measurements which mi ght be an indication 
that t he special emissivity model for Vesta (Mi iller & Lagerros! 
1998) needs a small adjustment for this wavelength range. 
It looks like t he emissivit y drop which has been seen by 
iRedman et all (fl99l [1998) occurs already at shorter wave- 



6 F. Hormuth and T. G. Miiller: Catalogue of ISO LWS observations of asteroids 

Table 4. Object specific TPM input parameters. The A ec i and /3 ec i values are the ecliptic coordinates of the spin vector direction, with 
j3 ec i counted from the equator. 



Object 


[H,G] 


Pm [h] 


'tec/, Peel 


source of shape & spin vector 


(1) Ceres 

(2) Pallas 
(4) Vesta 

(10) Hygiea 


[3.28,0.05]' 
[4.13,0.16] 2 
[3.20,0.34] 3 
[5.43,0.07] 2 


9.074184 
7.813225 
5.342129 
27.623265 


331.51°, +77.88° 
34.80°, -11.75° 
319.45°, +59.23° 
117.01°, -28.53° 


HST observations 4 

Lightcurve inversion techniques 5 - 6 

HST observations 7 

Lightcurve inversion techniques 8 



' lLagerkvist etal! dl992h; 2 iLaserkyist et all d200lh. . 
22/Sep/2005: 6 lTorppa et alj d2003h : 7 lThomas et alj lll997) 



Miiller & Lagerros (1998); 4 iThomas et al] d2005h : 5 M. Kaasalainen, priv. comm. 



M. Kaasalainen, priv. comm. 31/Jan/2006; 



lengths below 50/mi. Almost all LWS observations of Vesta 
were taken on one day (revolution 805) under very similar as- 
pect angles. This might also influence the outcome of the TPM 
technique which works best when combining data from differ- 
ent wavelengths, phase angles r otational phases and aspect an- 
gles dMiiller & Lagerrosll2002l) . IShevchenko & Tedescol {2006) 
derived an albedo of p# =0.370 from occultation measurements, 
while we obtained a radiometric value of py = 0.32+0.02, again 
this discrepancy might be explained by emissivity effects or the 
limited aspect angle range. The albedo value in Tb l. [5] agrees 
very well with previous studies by Miiller & Lagerros ( 1998J) (pv 
= 0.33) based on a set of thermal observations from mid-IR to 
mm- wavelengths . 

Our result for (10) Hygiea deviates significantly from 
previous studies: iTedesco et al.l d2002bl) found values of 
407. 12 +6.8 km and p v =0.0717+0.002. iMiiller & Lagerros! 
d 1998b gave 429.9 km and 0.066. But Hygiea's shape model 
is not that well defined due to a poor lightcurve coverage 
(in rotational phases) and some very low quality lightcurves 
(Kaasalainen, priv. comm.). There might also be a second pole 
solution at around A ec i = 300° (instead of the 117° used here). 
But neither using the second pole solutions nor using a spherical 
shape model lower the predicted size values (or the standard 
deviations) significantly. One possibility for the discrepancies 
in the diameters could be emissivity issues. If Hygiea's far-IR 
emissivity is lower then the corresponding diameter would also 
be smaller (and the albedo a bit larger). But the data quality 
and the very limited aspect angle range (all LWS observations 
were taken within 20 days) are not sufficient to draw firm con- 
clusions. Especi ally when looking at th e best quality direct size 
measurement by Ragazzoni et al. (2000): They gave an effective 
diameter of 444+35km (with an axis-ratio of a/b=l.ll), based 
on speckle techniques. Our radiometric results are still well 
within these error bars. 



5.3. Far-IR thermal lightcurve of (4) Vesta in the range 
45-1 05 fxm 

The visual lightcurve of (4) V esta is dominated by the influ- 
ence of the albedo variations (Dege wii et al.l [19791) . Standard 
lightcurve inversion techniques failed to produce a reliable 
shape model (Kaasalainen, priv. comm.), but high resolution 
HST imaging allowed a solution for the s hape and spin vector 
dThomas et al.lll997l) . lRedman et al.l dl992l) found that the 1 mm 
light-curve is apparently dominated by the triaxial shape, with- 
out any significant cont ributions from the optical albedo spots. 
Miiller & Barnes (2007) showed that at even longer wavelengths 
(around 3.2 mm) the mm-lightcurve follows for a large fraction 
of the rotational period the shape-introduced variations. They 
also demonstrated that the rotational phases with clear devia- 



tions are connected to structures (e.g., the Olbers feature) whi ch 
are visible in the HST images of (4) Vesta dBinzel et alJI 19971) . 

Before we combined TPM lightcurve predictions with the 
observed fluxes, we explored the influence of TPM input param- 
eters on the thermal lightcurve amplitude and phase for the given 
wavelength range and aspect angles. It turned out that neither 
the thermal inertia (the range between and 50 J irT 2 s~ 05 K -1 
was considered) nor the beaming parameter, describing various 
roughness scenarios, influenced the predictions for our case sig- 
nificantly. The predicted lightcurve amplitudes varied between 
5-6% (peak-to-peak) and the lightcurve phases only by a few 
minutes between the low and high thermal inertia predictions. 
Our data sets were therefore not sufficient to confine any of these 
parameters. 

The Vesta L02 observations from revolution 805 (28-Jan- 
1998) were taken over a period of about 6 hours, covering 
slightly more than one full rotation period of 5.34 hours. On the 
basis of the shape and spin-vector solution given in Tbl. [4] and 
the diameter and albedo results from Tbl. [5] we predicted the 
thermal lightcurves at the L02 key-wavelengths for the LWS ob- 
servations (note that the reference time frame is that of the ISO 
satellite and not the asteroid-intrinsic one). 

Figure [3] shows the observations at 102.3/mi together with 
the model prediction on an absolute flux and time scale. Figure|4] 
shows the combined data below 105 fim, first normalised per de- 
tector over the full observing period and then averaged for a 
given rotational phase angle. The quality of the data from the 
longer wavelength channels beyond 105yum was not sufficient 
for this kind of analysis. 

Both figures show that the observations follow the model 
predictions. The errorbars in Fig. |4] correspond to the error of 
the mean values plus a contribution from the flux uncertainty 
in these detectors. This figure combines in total 78 indepen- 
dent measurements (6 detectors x 1 3 measurements at different 
rotational phases). The model predictions match the observed 
lightcurve amplitude as well as the lightcurve phase (in the 
observer's time frame). We also indicated the rotational phase 
where the Olbers-feature (IZellner et al.lll997blBinzel et al.ll 19971) 
on Vesta's surface would be visible. We can not see any influ- 
ence o f this surface feature with lower albedo. IMiiller & Barnesl 
(|2007|) presented a more detailed discussion on how various ther- 
mal properties of this surface feature would influence the thermal 
lightcurve variations at different wavelengths. They found strong 
emissivity variations which they attributed to the Olbers struc- 
ture and a neighbouring region. The deviation from the model 
predictions at rotational phases between 240-280° in Fig.|4]are 
also close to the Olbers region and might therefore be the far-IR 
signature of the ejecta material which is deposi ted on one side 
of the Olbers impact structure as discussed in IMiiller & Barnesl 
(120071) . 
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Vesta: L02-LW1 obs. and TPM 



Q 




1 2 3 4 5 6 

Time [hours] since 28-Jan-1998 10:00 UT 

Fig. 3. The observed (dots with error bars) and the predicted ther- 
mal lightcurve (continuous line) of (4) Vesta at 102.3 jum (detec- 
tor LW1). 



Vesta: L02-SW1...LW1 obs. and TPM 



1.10 



1.05 



1.00 



0.95 



0.90 




100 200 
Rotational phase [deg] 



300 



Fig. 4. The observed (dots with error bars, normalised and av- 
eraged over five detectors (SW1, SW2, SW3, SW4, SW5 and 
LW1) and the predicted thermal lightcurve (continuous line) of 
(4) Vesta as a function of the rotational phase. 



6. Conclusions 

Via radiometric techniques we derived the first far-IR-based di- 
ameter and albedo values for the four large asteroids Ceres, 
Pallas, Vesta, and Hygiea. Based on photometry in the wave- 
length range 46-85/mi we calculated effective diameters D e ff 
of 960+17km, 534+15km, 549±13km, and 469±27km for 
these four objects. The corresponding albedo values py are 
0.096±0.003, 0.142±0.009, 0.317+0.015, and 0.056+0.006, re- 
spectively. We found excellent agreement with direct measure- 
ments from HST (Ceres and Vesta) and speckle/occultation re- 
sults (Pallas and Hygiea). Our data contain also the first far- 
IR lightcurve of an asteroid: Vesta's thermal lightcurve follows 



nicely the shape-introduced variations. Albedo and grain effects 
are important in the visual and mm-lightcurves, but seem to play 
a very small role at far-IR wavelengths. 

The catalogue of LWS observations of asteroids contains 
photometric and spectroscopic measurements of the main repre- 
sentatives of the main asteroid belt. The data will remain useful 
for thermophysical and mineralogical studies of these bodies, es- 
pecially in the context of the Dawn space mission to Ceres and 
Vesta. 
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Table 1. Overview of standard LWS asteroid observations in our catalogue. For the on-source pointings the distance to the asteroid A, the Sun- 
asteroid distance r, and the phase angle a are given, based on MPC (Minor Planet Center) ephemerides. The TDT number is a unique identifer, 
allowing to find the observation in the archive. 



Object 




II vi 

1D1 


Mid-observation UTC 


A(J1 


A [AU] 


r[AU] 


a [°] Remarks 


(1) 


Ceres 


09300401 


iyyo-uz-io 


U4. jl.lo 


L02 


2.673 


2.695 


21.2 






09304102 


iyyo-uz-i o 


1 J. 1 j.Zj 


L02 


(off-source reference measurement for TDT 09300401) 






10500402 


iyyo-uj -ui 


Uj.ZZ.ZZ 


L02 


2.520 


2.705 


21.5 






11900214 


iyyo-uo-1 j 


vZAZ. j 1 


L02 


2.344 


2.715 


21.1 






11905611 


iyyo-uj-i j 


10. j4.40 


L02 


(off-source reference measurement for TDT 1 1900214) 






12600114 




UZ.Uj.j4 


L02 


2.258 


2.721 


20.5 






25800302 


iyyo-u/-j i 


ly: 30:2/ 


L02 


2.258 


2.823 


19.2 






25805903 


iyyo-uo -ui 


Uo. jj.j4 


L02 


(off-source reference measurement for TDT 25800302) 






26500301 


iyyo-ua-u / 


in,] r\.'*>'7 
ly. 1U.Z/ 


L02 


2.351 


2.828 


20.0 






26505602 


iyyo-Uo -Uo 


uy: jZ:jo 


L02 


(off-source reference measurement for TDT 26500301) 






32100204 


iyyo-iu-uz 


1 J. Jl.JJ 


L02 


3.126 


2.868 


18.6 






32103506 


iyyo-iu-uj 


Uj.jZ.44 


L02 


(off-source reference measurement for TDT 32100204) 






53802209 




1 7-1 fvOQ 

1 Z,. 1U.U7 


L02 


3.112 


2.971 


18.9 






57902409 


1997-06-17 


09:56:26 


L02 


2.575 


2.978 


19.4 






59401908 


1997-07-02 


06:42:05 


L02 


2.391 


2.980 


17.9 






72001901 


1997-11-05 


01:10:04 


L02 


2.522 


2.973 


18.6 






74803304 


1997-12-03 


03:04:58 


L02 


2.901 


2.966 


19.3 






74803403 


1997-12-03 


08:11:59 


L01 


2.904 


2.966 


19.3 contains warmup artefacts 






75502902 


1997-12-10 


00:30:08 


L01 


(off-source reference measurement for TDTs 74803...) 






75503003 


1997-12-10 


01:38:00 


L02 


2.994 


2.964 


19.0 






76200502 


1997-12-16 


13:29:15 


L02 


3.081 


2.962 


18.6 






76903102 


1997-12-24 


01:08:50 


L02 


3.177 


2.960 


18.0 






76903203 


1997-12-24 


02:17:12 


L01 


3.178 


2.960 


18.0 


(2) 


Pallas 


23000306 


1996-07-03 


21:24:01 


L02 


2.497 


2.823 


20.9 






23002907 


1996-07-04 


12:01:22 


L02 


(off-source reference measurement for TDT 23000306) 






25100202 


1996-07-24 


19:49:12 


L02 


2.780 


2.871 


20.6 






25103603 


1996-07-25 


09:58:26 


L02 


(off-source reference measurement for TDT 25100202) 






26500503 


1996-08-07 


19:43:17 


L02 


2.967 


2.903 


19.9 






26505204 


1996-08-08 


07:27:57 


L02 


(off-source reference measurement for TDT 26500503) 






27200203 


1996-08-14 


18:19:41 


L02 


3.057 


2.919 


19.3 






27202004 


1996-08-15 


03:27:01 


L02 


(off-source reference measurement for TDT 27200203) 


(4) 


Vesta 


24402202 


iyyo-u/-i / 




L02 


1.593 


2.149 


26.6 






24404603 


iyyo-u/-io 


Uo:zy:4y 


L02 


(off-source reference measurement for TDT 24402202) 






80500101 


1 QQQ f\1 OQ 

iyyo-ui-zo 


lU.JO.Zj 


L02 


2.566 


2.549 


22.2 






80500104 


1 QQQ f\1 OQ 

iyyo-ui-zo 


1 1 .UZ.jU 


L02 


2.566 


2.549 


22.2 






80500107 


1 998-01 -28 

1770 V/l i-O 


1 1:29:15 


L02 


2.567 


2.549 


22.2 






80500110 


1998-01-28 


11:55:40 


L02 


2.567 


2.549 


22.2 






80500113 


1998-01-28 


12:22:05 


L02 


2.567 


2.549 


22.2 






80500116 


1998-01-28 


12:48:30 


L02 


2.567 


2.549 


22.2 






80500119 


1998-01-28 


13:14:55 


L02 


2.568 


2.549 


22.2 






80500122 


1998-01-28 


13:41:20 


L02 


2.568 


2.549 


22.2 






80500125 


1998-01-28 


14:07:45 


L02 


2.568 


2.549 


22.2 






80500128 


1998-01-28 


14:34:10 


L02 


2.568 


2.549 


22.2 






80500131 


1998-01-28 


15:00:35 


L02 


2.569 


2.549 


22.2 






80500134 


1998-01-28 


15:27:00 


L02 


2.569 


2.549 


22.2 






80500137 


1998-01-28 


15:53:25 


L02 


2.569 


2.549 


22.2 


(10) 


Hygiea 


83201702 


1998-02-24 


22:08:03 


L01 


3.179 


3.444 


16.6 CII background emission 






83201803 


1998-02-24 


22:44:05 


L02 


3.179 


3.444 


16.6 






84801302 


1998-03-12 


16:13:06 


L02 


3.402 


3.434 


16.7 






85303402 


1998-03-17 


16:09:46 


L02 


3.472 


3.430 


16.6 
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Table 2. Overview of L99 asteroid observations in our catalogue. All observations are LOl-like scans covering the full wavelength range. 



Object 


TDT 


Mid-observation UTC 


A [AU] 


r [AU] 


a [°] 


Remarks 


(1) Ceres 


07500601 


1996-01-31 


07:52:11 


2.891 


2.682 


19.9 


/—ITT 11 1 * 

CII background emission 




07500701 


1996-01-31 


10:10:56 


2.890 


2.682 


19.9 


CII background emission 




07501301 


1996-01-31 


16:08:13 


2.887 


2.682 


19.9 


CII background emission, 
















non-standard bias voltages 




07501401 


1996-01-31 


18:32:03 


2.886 


2.682 


20.0 


CII background emission, 
















non-standard bias voltages 


(4) Vesta 


05500201 


1996-01-11 


07:10:35 


2.313 


2.250 


24.8 


non-standard bias voltages 




05500301 


1996-01-11 


08:54:29 


2.312 


2.250 


24.9 


non-standard bias voltages 




05500401 


1996-01-11 


10:39:34 


2.311 


2.250 


24.9 


non-standard bias voltages 
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List of Objects 

'(1) Ceres' onpageQ] 
'(2) Pallas' on pagejl] 
'(4) Vesta' on page [I] 
'(10) Hygiea' on page[l] 
'Uranus' on page|j 



